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Abstract 


The  properties  of  the  impedance  and  scattering  matrix 
describing  waveguide  discontinuities  are  examined;  both  propaga¬ 
ting  and  evanescent  modes  are  considei*ed. 

It  is  shown  how  different  normalization  conditions  for  the 
normal  mode  solutions  in  the  guide  affect  the  impedance  matrix, 

A  suitable  choice  of  normalization  always  leads  to  a  symmetric 
imaginary  impedance  matrix  for  a  lossless  structure. 

The  scattering  matrix  is  no  longer  symmetric  or  unitary.  The 
simple  relationship  S  -  (Z  -  U)(Z  ♦  U)  1  is  shown  to  hold  only 
under  special  normalization  conditions. 

Next  the  matrices  describing  a  plane  of  lossless  obstacles 
arranged  in  a  periodic  array  are  examined.  A  different  type  of 
normalization  condition  must  be  used  here,  since  the  normal  modes 
are  orthogonal  in  the  conjugate  sense  (bi-orthogonal) . 

Although  the  structure  is  reciprocal,  none  of  the  matrices 
are  symmetric.  A  suitable  normalization  leads  to  a  skev-hermitian 
impedance  matrix  and  to  a  unitary  submatrix  of  the  scattering 
matrix  corresponding  to  propagating  modes. 


Matrix  Description  of  Waveguide  Discontinuities  l 
Plane  Arrays  of  Scatterers  in  the  Presence  of  Evanescent  Modes 


In  treating  the  problem  of  propagation  of  electromagnetic 
waves  past  waveguide  discontinuities  it  is  often  convenient  to 
define  a  set  of  equivalent  voltages  and  currents,  corresponding  to 
linear  combinations  of  incident  and  reflected  wave  amplitudes,  there¬ 
upon  reducing  it  to  a  circuit  problem.  The  equivalent  circuit  fcr 
the  discontinuity  lends  itself  to  description  in  terms  of  the  usual 
circuit  type  matrices  such  as  the  impedance,  admittance,  scattering 
or  other  matrix. 

The  properties  of  the  above  matrices  have  been  described  in 
1  2 

the  literature  9  when  only  propagating  modes  are  considered.  This 
is  normally  the  case  when  the  terminal  planes  on  which  voltages  and 
currents  are  defined  are  chosen  far  enough  from  the  discontinuity. 
However  in  many  cases  one  cannot  neglect  the  effect  of  the  evanescent 
modes,  as  for  example  when  two  discontinuities  are  closely  spaced. 

The  purpose  of  this  report  is  to  study  the  properties  of  the 
matrices  describing  waveguide  discontinuities,  or  plane  lattices 


1.  C.G.  Montgomery,  R.H.  Dicke  and  E.M.  Purcell;  "Principles 
of  Microwave  Circuits",  vol.  8  of  MIT  Rad.  Lab.  Series, 
McGraw-Hill  Book  Company,  Inc.,  New  York,  1948. 

2.  D.  Kerns:  "Basis  of  Application  of  Network  Equations  to 
Waveguide  Problems",  Journal  of  Research  NBS,  pp.  515-540, 
May  1949, 
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of  scatterers,  when  both  propagating  and  evanescent  modes  must  be 
considered. 

It  is  important  to  note  that  voltages  and  currents  are 
defined  quantities  and  as  such  they  can  be  chosen  in  different 
ways.  Correspondingly  the  properties  of  the  resulting  matrices 
will  be  affected  by  that  choice.  An  alternate  way  to  see  this 
arbitrariness  in  the  definition  of  the  voltages  and  currents  is  to 
note  that  the  normal  mode  solutions  in  the  waveguide  can  be  nor¬ 
malized  in  different  ways. 

To  be  specific  consider  a  rectangular  waveguide  with  a  dis¬ 
continuity  extending  in  the  Mz"  direction  from  z^  to  z^  as  in 
Fig.  1.  The  transverse  fields  on  the  two  sides  of  the  discontinuity 


i,  i 


I 

I 


Fjg«  1  A  general  waveguide  discontinuity 
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can  be  expanded  in  a  set  of  normal  modes: 
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for  z  <  z„ 


r»<,-*2) 


e  and  h  are  normal  mode  functions  in  the  guide  and  are  related 
n  n 

by  means  of  a  dyadic  impedance  or  admittance: 


where 
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(2b) 
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R.E.  Collin:  "Field  Theory  of  Guided  Waves",  Chapter  5, 
McGraw-Hill  Book  Company,  Inc.,  New  York,  1960. 
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Z  and  Y  arc  the  scalar  wave  impedance  and  admittance  and  are 
n  n 

real  for  propagating  modes  and  imaginary  for  evanescent  modes. 

^  a 

The  normal  mode  functions  e^,  h^  are  orthogonal  in  a  waveguide. 
Also  one  of  the  two  functions  can  always  be  chosen  real.  We  will 
assume  in  the  following  that  e^  is  real.  These  functions  can  be 
normalized  in  a  variety  of  ways,  some  of  the  possibilities  will 
now  be  considered. 

Let : 


e  dS 
m 


Y  N  6 
m  n  nm 


(3) 


where  S  denotes  the  waveguide  cross  section  and  is  a  real  posi¬ 
tive  normalization  constant,  arbitrary  as  yet. 

We  choose  the  voltages  to  be  proportional  to  the  amplitude 
of  the  transverse  electric  field  and  the  currents  proportional  to 
the  amplitude  of  the  transverse  magnetic  field,  thus: 


/  -  K.  a  V"  -  Lb 

n  in  n  n  In  n 


(4) 
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In  order  to  keep  the  complex  power  flow  invariant  it  is  necessary 


that : 
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(5b) 


Also  we  may  choose: 


v 

n  m  In 


n 


2n 


cn 


(6) 


where  Z  is  any  convenient  characteristic  impedance  for  the  equiva¬ 
lent  transmission  line. 

It  can  be  seen  from  (5b)  that  for  propagating  modes  K^n  and 

K~  can  both  be  chosen  real  since  Y  is  real.  However,  for  evanes- 
<n  n 

cent  modes  at  least  one  of  the  two  contacts  must  be  chosen  imaginary. 

One  can  still  choose  both  K^n  and  real  in  all  cases  pro¬ 

vided  the  definition  of  complex  power  flow  is  modified  as  follows: 


I  vO*  ■  I  y„i*/  ■  I  vLi*/ 


(7a) 


for  propagating  modes 


i  I JVO*  ‘  \  '  i  \  <7b> 


for  evanescent  modes 
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I n  (7b)  the  upper  sign  holds  for  H  inodes  and  the  lower  sign 
for  E  modes. 

Before  actually  choosing  a  specific  normalization  we  will 
first  derive  some  general  properties  of  the  impedance  matrix. 

We  will  number  the  modes  on  the  two  sides  of  the  discontinuity 
in  consecutive  order,  that  is,  we  can  define  a  voltage  vector  and 
current  vector  (column  matrix): 


v  -  v*  +  v- 

i  -  i*  -  r 


(8) 


where 


V~  - 


and  similarly  for  the  currents.  In  view  of  the  linearity  of 
Maxwell's  equations  we  have: 


V  -  ZI 


(9) 


where  Z  is  the  impedance  matrix. 

Consider  first  two  independent  solutions  to  Maxwell's  equations 
satisfying  boundary  conditions  in  the  guide.  The  following  relation 
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then  holds  in  a  reciprocal  medium. 


(E1 


x  H2 


E2  x  H1) 


(10) 


the  superscripts  refer  to  the  two  independent  solutions.  When  this 
relation  is  integrated  over  a  region  enclosing  the  discontinuity 
we  obtain: 


(I1  ft1)  dv 


,  p  9  ,  Y  N 
(V1!2  -  V2!1)  n  n 
n  n  n  n 


2n 


r  r  x  ft2  -  ft2  x  ft1) .  dft 


(id 


where  it  was  assumed  that  the  transverse  fields  have  been  expanded 
in  normal  modes  with  the  normalization  as  given  by  (3). 

Consider  now  two  conditions  for  the  terminal  planes: 


Condition.  1.  V1  -  0  n  /  1  I2  -  Y^V2 


Condition  2,  V2  -  0  n  /  j  I1  -  Y^V^ 


Then  it  follows  from  (11): 


.  Vi 

iJ  ^li^  Ji  ^lj*2J 


(12) 
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where  Y  and  Y  are  elements  of  the  admittance  matrix.  The 
admittance  matrix,  and  thus  the  impedance  matrix,  will  be  symmetric 
if: 


Vi 


2i 


Y1N.1 

*1/23 


(13) 


At  this  point  we  will  make  the  choice  of  normalization. 
Consider  the  following  cases: 

Case  It  f  f  e _  •  e_  -  N_  -  1  (14) 

-  J  £  n  n  n 

Using  condition  (5b)  we  must  have: 


hn'ln 


Y* 

n 


We  can  then  choose  ■  1  and  “  ^n.  With  this  choice 
relation  (13)  is  satisfied  and  thus: 


Case  II: 


e  dS 

n 


and  also 


for  propagating  modes 
for  evanescent  inodes 


(IS) 


n 

The  (+)  and  ( -)  signs  correspond  to  E  and  H  modes.  Again  using 
condition  (5b)  we  must  have: 
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^ln^2n  ”  ^  for  Pr0Pa8*lting  modes 

#  — 

^ln^2n  "  f°r  evanescen^  rood® s 

This  then  allows  us  to  choose: 

^ln  "  *gn  "  1  for  propagating  modes 

K^n  *  l^K^  "  — 1 J  for  ®van®8cent  modes 

With  that  choice  relation  (13)  is  satisfied  and  both  the 
admittance  and  impedance  matrix  are  again  symmetric* 

Further  properties  of  the  impedance  matrix  can  be  derived 
from  the  energy  condition.  Integrating  the  complex  Poynting 
vector  over  a  region  containing  the  discontinuity  we  obtain: 


P  ♦  2jo)(W  -  W  ) 

m  e 


where  the  tilda  stands  for  the  transposed  matrix.  If  the  structure 
is  lossless: 


P  -  0  -  |  Re (VI*)  -  i  £l*(Z  ♦  Z*)  I] 


Z  ♦  Z*  -  0 


Then 


(17) 
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sin  ce  Z  •  Z. 

We  conclude  that  under  normalization  condition  (14)  or  (15) 
with  the  definition  of  complex  power  in  t he  sense  of  (5b)  the 
impedance  matrix  of  a  lossless  discontinuity  is  symmetric  and 
imaginary.  The  same  holds  for  the  admittance  matrix. 

It  is  sometimes  convenient  to  choose  the  characteristic 
impedances  of  the  equivalent  transmission  lines  as  unity  for 
all  modes.  As  we  have  seen  before  condition  (5b)  does  not  allow 
us  that  choice  unless  we  resort  to  the  modified  definition  of 
complex  power  flow  as  given  by  (7a)  and  (7b).  Using  normalization 
condition  (15)  we  obtain: 


i 

i 


- 1 


for  both  propagating  and  evanescent  modes.  It  is  now  possible 
to  choose 


K 


In 


-  1 


The  properties  of  the  impedance  matrix  are  now  somewhat 
different.  It  is  convenient  to  partition  the  voltage,  current 
and  impedance  matrices  as  follows: 


- 1 

< 

_ 1 

- 1 
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vhere  the  subscripts  p  and  e  refer  to  propagating  and  evanescent 
modes. 

The  reciprocity  relation  (11)  now  becomes: 

V1  I2  -  V2  I1  ♦  jr?  i2  -  V2  I1!  -  0  (19) 

PP  p  p  -  u  e  e  e  e-* 

The  upper  sign  holds  for  E  inodes  and  the  lower  sign  for  H  modes. 
Consider  now  the  following  conditions  for  the  terminal  plane sj 

Condition  1  T  I1!- 

-  L  p-l 


With: 


Ji1-  P 


we  obtain: 


The  submatrix  £z  ]  is  symmetric.  We  now  assume  that 
only  E  or  H  modes  are  present. 


Condition  2 


we  then  obtain 


Condition  3 


With: 


we  obtain: 


Z_  -  ZJ,  for  the  V  2  ”1  submatrix. 

ij  ji  L  ee-1 


■  0 

pe  -  J  ep 


or  more  generally: 
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If  we  now  impose  the  energy  condition  for  a  lossless 
structure  on  the  impedance  matrix  we  have: 


P  -  0  -  Ref  I*v  ♦  J?V  1  «  Re!"  I*z  I  ] 

L  P  P  ~  «  e-1  L  PPPPJ 


♦  Re[>jI*Z  I  1  ♦  Refl*Z  I  1  ♦  Ref>jI*Z  I  H 
u  -  e  ee  e-J  *-  p  pe  eJ  1 — °  e  ep  p-* 


Therefore : 


C  ZppD  ♦  C  ZppD  -  0  C  zppD  imaginary 


CZeeD  -  EZeeD  -0  EZeelisreal 


Also : 


rr  z  ♦  jz*  ii  » o 

p*-  pe  J  ep->  e 


Therefore : 


but : 


K 7  0 

&p. :  ° 


Then : 


r  z*  *  z  i-o 

*-  ep  ep-1 


(21) 
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Therefore 


£  Z  J  is  imaginary 

I"Z  1  is  real 
*-  pe-1 


An  alternate  description  of  waveguide  discontinuities  is  by 
means  of  the  scattering  matrix. 

We  define  the  matrix  by: 


V"  -  SV*  (22) 

If  we  choose  1C  ■  Kn  *1  for  all  modes  then: 
in  2n 

v  -  V*  ♦  v'  -  Z(I*  ♦  l“)  -  Z(V*  -  V*) 

(Z  -  U)V4  -  (Z  ♦  u)v“ 

and  S  -  (Z  ♦  U)'1  (Z  -  U)  (23) 

The  matrix  S  given  by  (23)  is  not  symmetric  since  Z  is  not 
symmetric  for  this  case.  Also  S  is  not  an  unitary  matrix  any  more. 

Condition  (14)  and  (15)  do  not  lead  to  a  simple  relationship 
between  the  impedance  matrix  and  the  scattering  matrix,  because 
the  characteristic  impedances  of  the  equivalent  transmission  lines 
are  real  for  propagating  modes  and  imaginary  for  evanescent  modes. 
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The  properties  of  the  scattering  matrix  as  given  by  (23)  will  be 
studied  in  a  later  part  of  this  report. 

We  now  consider  a  plane  of  lossless  obstacles  arranged  in  a 
periodic  array j  no  longer  is  it  possible  here  to  use  normalization 
conditions  (14)  or  (15)  because  the  normal  mode  solutions  in  the 
structure  are  always  complex  such  as  to  satisfy  periodic  boundary 
conditions  in  the  transverse  plane. 

The  following  normalization  procedure,  in  the  conjugate  sense 
can  be  used  for  both  a  lossless  discontinuity  in  a  guide  and  an 
array  of  lossless  scatterers. 

Let : 


m 


e  dS 
m 


Y*N  6 


m  n  run 


(24) 


where  S  denotes  the  cross  section  of  the  unit  cell  in  the  structure 

and  N  is  a  oositive  real  normalization  constant, 
n 

It  should  be  noted  that  it  is  still  possible  to  normalize 
the  modes  in  a  lossless  periodic  structure  as  in  (14)  or  (15) 
provided  one  restricts  propagation  to  incidence  normal  to  the 
plane.  That  is  indeed  the  case  since  the  E  and  H  modes  are 
orthogonal  in  a  periodic  structure  of  this  sort  just  like  in  a  wave 
guide.  Again  consider  two  cases: 


p  r  _ 

e  •  e  dS 
n  n 

b 


JJ 


Case  I 


(25) 
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Case  II 


I 


n  ^ 


•  e  dS 

n 


for  propagating  modes 
for  evanescent  modes 


(26) 


The  upper  sign  holds  for  E  modes  and  the  lower  for  H  modes.  A 
reciprocity  relation  can  be  written  for  the  field  in  the  struc¬ 
ture  as  follows : 


r>  •  (E1  x  H*2  ♦  £*2  xt1)  -  0  (27) 

However,  this  relation  assumes  a  lossless  structure.  There¬ 
fore,  it  is  easier  to  use  directly  the  energy  condition. 

The  complex  Poynting  vector  integrated  over  the  terminal 
planes  of  a  unit  cell  yields 

|  f  t  x  ft*  ♦  <1$  -  2j0)(Wm  -  We)  -  |  VI*  (28) 

if  one  defines  complex  power  in  the  sense  of  (5b)  with  normalization 
given  by  (25)  or  (26). 

The  structure  is  lossless  so 

P  -  0  »  -  He (VI  )  -  i  (I  ZI  ♦  I  Z  I) 
or  Z  ♦  Z  ■  0.  The  impedance  matrix  is  thus  skew  hermitian. 
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However,  if  normalization  condition  (26)  is  used  together 
with  the  power  definition  (7a)  and  (7b)  the  complex  integration 
yields : 


Re(I*Z  I  )  ♦  Re(*JI*Z  I  )  ♦  Re(I*Z  I  )  ♦  fte(+j?Z  l)  -  0 
p  pp  p  — J  e  ee  e  t>  ds  e  —  e  eD  d 


p  pe  e' 


e  ep  p 


where  £  is  partitioned  as  in  (18).  We  moat  then  have: 


or  |Z  I  is  skew  -  hermitian 

LppJ 


or  Z  is  hermitian 

Le«J 


This  impedance  matrix  is  related  to  the  scattering  matrix  by 
means  of  relation  (23).  No  simple  relationship  exists  between  the 
S  and  Z  matrices  obtained  from  any  other  normalization. 

To  obtain  the  properties  of  the  scattering  matrix  we  partition 
it  as  follows: 


PH 
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'M 
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_  e_ 

J>] 
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V® 
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(29) 


The  real  power  is  given  by: 
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P  -  0 


<*;>*  »;  -  v  v  i  *;  -  <v*  o 


-  0 

(30) 


for  a  lossless  structure. 
Using  (29)  this  becomes: 


~  * 


(V*)  ru  -  s*  s  iv*  -  ( v *)  rs*  s 

'  P  L  PP  pp-1  P  e'  *-  pe  pe 


♦  J(S 
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e'  L  pe  pp 


IS  Iv*  -  (v*)*  rs*  s  ♦  js*  Iv*  «  o 

J  epJ  p  v  p'  *-  pp  pe  -  J  ep-1  e 


This  then  yields: 


°r  y 18  — 

[s«  -  5J*  i>RW 
EpS*]*Os*rI"  0 

The  above  scattering  matrix  is  defined  with  normalization  given  by 
(26)  and  power  definition  given  by  (7a)  and  (7b) . 

In  conclusion  we  have  shown  that  for  a  lossless  discontinuity 
in  a  waveguide  it  is  always  possible  to  choose  the  equivalent 
voltages  and  currents  such  as  to  obtain  an  impedance  matrix  which 
is  symmetric  and  imaginary  just  like  in  the  case  when  only  propaga¬ 
ting  modes  are  present.  For  this  case  however  the  characteristic 


impedances  of  the  equivalent  transmission  lines  are  real  for  propa¬ 
gating  modes  and  imaginary  for  the  evanescent  modes. 

If  one  tries  to  make  all  characteristic  impedances  of  the 
transmission  lines  real  for  all  modes  the  corresponding  impedance 
matrix  is  no  longer  symmetric  or  imaginary. 

For  a  periodic  array  of  scatterers  it  is  possible  to  choose 
voltages  and  currents  such  that  the  impedance  matrix  is  skew-hermit ian. 
Although  the  structure  is  reciprocal  the  impedance  matrix  is  not 
symmetric. 

The  scattering  matrix  is  no  longer  symmetric  or  unitary  in 
either  case;  it  is  only  possible  to  obtain  a  unitary  subtnatrix 
corresponding  to  propagating  modes. 


